Atom-Pair Tunneling and Quantum Phase Transition in Strong 

Interaction Regime 

J.-Q. LiangQ J.-L. Liu, W.-D. Li, and Z.-J. L 

Institute of Theoretical Physics and Department of Physics, 
Shanxi University, Taiyuan, 030006, China 

We propose a Hamiltonian of ultracold spinless atoms in optical lattices including 
the two-body interaction of nearest neighbors, which reduces to the Bose-Hubbard 
model in weak interaction limit. An atom-pair hoping term appearing in the new 
Hamiltonian explains naturally the recent experimental observation of correlated 
tunneling in a double-well trap with strong atom-atom interactions and moreover 
leads to a new dynamic process of atom-pair tunneling where strongly interacting 
atoms can tunnel back and forth as a fragmented pair. Finally a new dynamics 
of oscillations induced by the atom-pair tunneling is found in the strong interac- 
tion regime, where the Bose-Hubbard model gives rise to the insulator state with 
fixed time-averaged value of atom-occupation-number only. Quantum phase transi- 
tions between two quantum phases characterized by degenerate and non-degenerate 
ground states are shown to be coinciding with the Landau second-order phase tran- 
sition theory. In the system of finite atom-number the degeneracy of ground states 
can be removed by quantum tunneling for the even-number of atoms but not for the 
odd-number. 

PACS numbers: 03.75.Lm, 37.10.Jk, 03.75.Kk, 75.75.+a 

I. INTRODUCTION 

The ultracold atom-gas clouds possess many advantages for investigation of quantum 
phenomena and hence become a test ground of quantum mechanical principles in many 
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extraordinary aspects and new prospective regime. Recently rapid advances of experimental 
techniques in optical traps open up a prospect for the study of quantum phase transition 
(QPT), for example, from a superffuid to a Mott-insulator[l, 2, 3], where the ratio between 
tunnel coupling through the interwell barriers and the atom-atom interaction plays a crucial 
role. Up to date the QPT has been studied only based on the well known and widely 
applied Bose-Hubbard (BH) Hamiltonian, which in a mean field description leads to coupled 
equations of population-imbalance and phase difference between two wells for a double- 
well trap and thus can describe the Josephson oscillation as well as self-trapping of Bose- 
Einstein condensates (BECs) qualitatively. It is demonstrated that the investigation of 
experimentally observed non-linear self-trapping (4] could provide a test ground of the mean 
field description of the BH model in the strong nonlinear regime. Particle tunneling through a 
classically prohibiting potential barrier is one of the characteristic effects of quantum physics, 
which has been essentially studied in a single-particle manner. For the many-body case 
strong interactions between particles may fundamentally alter the tunnel configuration and 
result in a correlated tunneling, which was explored most recently in ultracold atoms j^, E|. 
On the other hand a co-tunneling regime can be achieved in coupled mesoscopic quantum 

u n 

dots, where separate electrons only tunnel in a correlated way [3, |8J. We point out in this 
paper that the well known BH Hamiltonian based on the hard-core interaction is not able to 
describe the dynamics of atom-pair tunneling, which is the dominant dynamic effect in strong 
interaction regime, since the correlated atom-pair tunneling requires long range correlation 
of wave functions. The BH Hamiltonian, which is valid in a relatively weak interaction 
regime, should be extended in the strong interacting re gim e to include the superexchange 
interactions between atoms on neighboring lattice sites PJlol]. A peculiar atom-pair hopping 
term appearing in the new Hamiltonian explains very well the recently reported experimental 
observation of correlated tunneling and moreover leads to new dynamics in superstrong 
interaction regime, which has not yet been explored. 

The Hamiltonian of second-quantization beyond the on-site approximation is seen to be 
(see appendix for derivation) 

i i i 

+(U 1 + U 2 ) n i n i+i + ~Y ^{a\a\a i+ ia i+ i + a\ +l a\ +1 aia,i) (1) 
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where aj is the zth-site boson annihilation operator and rij = ajaj is the corresponding atom- 
number operator. The new coupling constants between atoms on neighboring lattice-sites 
are given by 

Ui= w i (x 2 )w i+ i(x 1 )U(\xi - x 2 \)wi +1 (xi)wi(x 2 )dx 1 dx2 1 



U 2 = / Wi(x 2 )wi(xi)U(\x 1 - x 2 \)w i+1 (xi)w i+1 (x 2 )dx x dx 2 , 

and 

U 3 = J Wi(x 2 )Wi + i(xi)U(\x 1 - x 2 \)w i+1 (x 1 )w i+ i(x 2 )dx 1 dx 2 , 

where U(\xi — x 2 \) is the two-atom interaction potential and Wi denotes Wannier wave 
functions. Uq is the usual on-site atom-interaction strength, and the fourth term relating to 
U\ denotes the nearest-neighbor repulsion, which has been given in various lattice models 
in the literature, while the U 2 part is new. The conceptually new transition-matrix element 
U 2 /2 describes obviously the atom-pair tunneling. The atom- atom interaction including 
the nearest-neighbor results in an additional Josephson tunneling term in relation with the 
coupling constant U 3 and the atom- number operator + 114+1 — 1. This term as a matter 
of fact suppresses the interwell hopping. The coupling constants Ui, U 2 , and U 3 are much 
smaller comparing with Uq and can be evaluated roughly with Gaussian Wannier function 
and 5-function potential. We have 



3 . 



where 



Ui « U 2 « e 2 U , U 3 ^eW , 



V 



with V , Er being the well depth and lattice recoil energy respectively. Using the experi- 
mental values^ for V and Er we find from the equations above that U\ ~ U 2 « 0.02f/ . 
This is in close agreement with the values we take later on with best fitting of the curves 
with experimental data points U\ — U 2 — 0.018t/ (see below). 

When a correlated many-body quantum-system is driven by a controllable parameter, 
the ground-state energy may have a structural change at a critical value of this parameter. 
This phenomenon is called QPT from a disordered phase to an ordered one obeying the 
Landau theory with ordinary symmetry breaking. A well-studied example for the QPT 
is a model described by the Hamiltonian of Lipkin-Meshkov-Glick (LMG) type, which is 
originally introduced in nuclear physics as exactly solvable many-body interacting model 
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and now has been used in a broad range of topics such as spin system and BECs. In the 
spin system QPT between ferromagnetic phase of long-range magnetic order and disordered 
paramagnetic phase can occur at a critical value of external magnetic filed. We show that 
the atom-pair tunneling term results in new dynamics of BEC in a double-well trap and 
a formally same QPT as that in the spin system can be found in the strong atom-atom 
interaction regime, where the BH Hamiltonian could only give rise to an insulator-phase 
with a fixed occupation number of atoms in each well. In this paper a new method is 
adopted for the investigation of QPT, where the many-body system is converted into a 
giant pseudo-spin which in turn is mapped to a single particle problem with an effective 
potential. Moreover for the finite-size system with a small number of particles the quantum 
tunneling may play a significant role and as a consequence the QPT would depend on the 
number-parity (even or odd) of particles resulted from the topological phase interference of 
tunneling paths. This peculiar phenomenon is analyzed in section IV. 



II. TWO-ATOM DYNAMICS IN A DOUBLE- WELL TRAP AND 
EXPERIMENTAL EVIDENCE OF ATOM-PAIR TUNNELING 

The tunneling dynamics of a few atoms loaded in a double-well trap has been studied 
by varying the interaction strength from weak to strong limit and it is shown for the two- 
atom case that the character of tunneling changes from Rabi oscillation to a correlated 
process with increasing interaction , namely, when repulsive interactions are strong, 
two atoms located on one side of the barrier cannot separate, but tunnel together as a 
pair in a co-tunneling process. As a matter of fact interactions of ultracold atoms can be 
adjusted experimentally over a wide range via Feshbach resonances which make it possible to 
explore the limit of strong correlation. A direct observation of the correlated tunneling was 
reported recently^] and theoretical analysis has been also presented in terms of two-body 
quantum mechanics (| . Moreover it has been demonstrated that if the bosons repel each 
other infinitely strongly, they can be mapped to noninteracting fermions (the fermionization 
limit) in the sense that the hard-core interaction mimics the exclusion principle. Near 
fermionization the strongly interacting atoms tunnel back and forth as a fragmented pair in 
a double- well trap jf|. 

For the double- well potential the Hamiltonian Eq.(l) reduces (with a constant-energy 
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term renormalized to zero) to 

H = y(^L - n R ) - J{a\a R + a R a L ) + ^ [n L (n L - 1) + n R (n R - 1)] 

+ (U 1 + U 2 )n L n R + -^-{a\a\a R a R + a R a R a L a L ), (2) 
where A = e L — e R is the bias potential between two wells and 

J = J - (N - 1)U 3 

with Jo being the single- atom Josephson coupling constant. The effective Josephson coupling 
constant J , which is suppressed by the nearest-neighbor repulsion coupling, can become a 
negative value for sufficiently large number of atoms N. In tight-confinement approximation, 
the quantum state of atoms can be described in a Fock state basis \ni,n R > with til and 
n R being non-negative integers. For one-atom occupation in the two wells corresponding 
to a filling factor 1/2 of the optical lattice, the Hamiltonian Eq.(2) gives rise to nothing 
new but the Josephson oscillations (see below). We are interested in the dynamics of two- 
atom occupation (filling factor 1 ), where the Fock-state base- vectors are |1, 1 >, |2, > and 
1 0, 2 >. The states |0,2 >, |2,0 > both couple to the state |1, 1 > by the Josephson-tunnel 
matrix-element J and couple each other directly via the matrix element U2, the effect of 
which is vanishingly small in weak interaction regime (roughly speaking Uq < J ), where 
the Josephson tunneling dominates. On the other hand for strong interactions (Uq > J) 
the energy difference between states |0,2 > and |1, 1 > is larger than the coupling J and 
therefore the Josephson-tunneling induced transition between these states is suppressed due 
to the strong detuning. However the direct transition between states |0,2 >, and |2,0 > 
via U2 matrix element is still resonant and gradually becomes the dominant dynamics with 
continuously increasing interactions. The BH Hamiltonian can only give rise to the first- 
order transitions between states |1,1 > and |0,2 > (or |2,0 >), while the transition, for 
example, from state |2, > to state |0, 2 > is generated only by a second-order transition^ 
via the J matrix element that |2,0 >^ |1,1 >^ |0,2 > with an effective coupling con- 
stant 2J 2 /Uq, which alone, we will see, is too small to be responsible for the correlated 
tunneling. We demonstrate in this paper that the correlated tunneling in fact includes both 
the second-order transition and the atom-pair tunneling via the matrix element U2, which 
is more important dynamic process in the superstrong interaction regime (Uo » J) near 
fermionization. To see this closely we can evaluate the energy eigenvalues (i — 1, 2, 3) and 
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the corresponding eigenstates \ipi > by direct diagonalization. The two-atom dynamics is 



studied by solving the time-dependent Schrddinger equation ih4:\ij)(t) >= H\tp(t) > within 
a three-state description for both the Eq.(2) and BH Hamiltonians. Thus the time-evolution 
of average position of the two atoms can be evaluated by 

< x >= - < ^{t)\n L -n R \ij){t) >, (3) 

which characterizes the tunneling dynamics of atoms. We consider the two atoms initially 
localized on one side of the double wells, for example the left-well. The time-evolutions 
of average positions are shown in Fig.l (solid lines) and are compared with the measured 
time-resolved traces (black dots). If the interaction energy Uq is much smaller than the 
tunnel coupling constant J (J /Uq = 1.5 in the experiment [5j), the Josephson tunneling is a 
dominant dynamic mechanism and the average-position oscillation contains more than one 
frequency component indicating the transition between the two states |2,0 > and |0,2 >, 
which is induced mainly by the second-order process of the Josephson-matrix element 
The direct transition between the two states |2, > and |0, 2 > via matrix element U2, which 
is more than two-order lower than the second-order Josephson tunneling, is negligibly small 
and thus results essentially in no effect seen from Fig. 1(a). However, when the interaction 
energy Uq increases to reach the interaction-dominated regime (J/Uq = 0.2 in the experiment 
[5I), where the Josephson tunneling is suppressed due to the energy difference between states 
|0, 2 > (or |2, >) and |1, 1 >, the main dynamic process is the correlated tunneling with fre- 
quency 550Hz (oscillation period 1.8ms) seen from the measured two-atom average position 
5| (Fig. 1(b) black dots), which coincides obviously with the numerical result of the Hamilto- 
nian Eq.(2) (red-solid line in Fig. 1(b)) including direct pair-tunneling process between states 
|0, 2 > and |2, > via the transition element U2, which leads to a significant modification of 
the dynamics. The interaction parameters U\ = U2 = 0.018[/o are determined by the best 

n 

fitting with the experimental data [5j for a fixed Josephson coupling constant Jo, which is 



obtained from the experiment of single-atom occupation [5j (see below). We see that the 
Josephson tunneling between the states |1, 1 > and |0,2 > (or |2,0 >) is visible as a small 
modulation with a period of 400/xs . The second-order transition |0,2 >— > |1, 1 >^ |2,0 > 
(or |2,0 >^ |1,1 >^ 1 , 2 >) with the coupling constant 2J 2 /Uq, which is more than 
four times of the U2 value in the present case, is indeed the dominant dynamic process as 
demonstrated in Ref . jsf . It was, however, realized in Ref . [5!] that the theoretical results eval- 
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uated from the BH Hamiltonian with the second-order transition give rise to slightly longer 
oscillation-period of average position and lower oscillation amplitude (Fig. 1(b), blue-dash 



line) than the experimental data (black dots [5]). This deviation although small is crucial and 
was compensated by modifying the Josephson coupling constant with an additional 3-10% 
higher value in Ref. |5j. We show that this problem can be cured by the atom-pair hoping 
via the transition element U%. When the interaction strength increases up to the value that 
J/Uq = 0.1, the atom-pair hoping becomes more important since the coupling constant of 
second-order transition is about the same value of U2 and thus the BH Hamiltonian fails to 
describe tunneling dynamics. The average positions of the two atoms evaluated from Eq. (2) 
(red-solid line) and BH Hamiltonian (blue-dash line) are given in Fig. 1(c) showing a great 
difference, that the oscillation frequency from Eq.(2) is almost two times higher than that 
evaluated from the BH Hamiltonian for the transition between states |0,2 > and |2,0 > (i. 
e. the oscillation of large amplitude in Fig. 1(c)). Increasing the atom-atom interaction the 
atom-pair hopping introduced in this paper becomes more and more important than the 
Josephson tunneling and at a critical point J = Jq — (N — l)Us = 0, where the Josephson 
tunneling vanishes, we have the sinusoidal population oscillations between two states |2, >, 
1 0, 2 > induced only by the atom-pair tunneling i.e. the effective Josephson oscillation of 
atom-pair. We see that atom-pair tunneling becomes the dominant dynamic process and 
should be observed experimentally. 

In addition to the average center-of mass positions the phase relation between two wells is 
also measured experimentally by a separate interferometric sequence of double-slit matter- 
wave interference pattern, which can be evaluated by 

p = J2 < a ia j > e ik < r '~ r ^ = N(l + Vcos[0 + k-(r L - r R )}) (4) 

where <fi is the relative phase between two wells, V is the visibility and iV denotes the 
total number of atoms in the double-well. For a single-atom occupation the experimentally 
observed signal is still sinusoidal however with a correspondingly lower frequency 5[] in the 
strong interaction regime (J/Uq = 0.2). Both Hamiltonians give rise to the same expected 
sinusoidal population oscillation between two states |1,0 >, |0, 1 > with oscillation-period 
depending only on the Josephson-matrix element J = Jq. As a comparison the corresponding 
average center-of mass position, phase and visibility are shown in Fig. 2 with black dots 
representing the experimental data [3] from which Josephson coupling constant J is obtained 
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Jo = 0.55 x 10 3 Hz- h, where h is the Planck's constant. We see that the two Hamiltonians 
make no difference in description of dynamic behavior of single-atom occupation seen from 
the coinciding red (solid) (proposed Hamiltonian) and blue (dash) (BH Hamiltonian) lines. 

III. ATOM- ATOM INTERACTION INDUCED DYNAMICS AND QUANTUM 



Dynamics of cold atoms in a double-well trap has been well studied based on the BH 
Hamiltonian. The 7r-phase oscillation, in which the time-averaged value of the phase differ- 
ence between two wells is equal to n, has been found in the weak interaction limit (J >> 
Uq). In the strong interaction region (Uq » J) the BH Hamiltonian results in the insulator- 
phase only, where the system is in the fixed atom occupation-number state. We show that 
the atom-pair tunneling in the Hamiltonian Eq.(2) can lead to a new dynamics of atom 
occupation-number oscillation, which has not yet been explored. We in this paper provide 
an analytic investigation based on an effective Hamiltonian of single-particle with canonical 
variables: the atom-number difference (or population imbalance) and phase difference be- 
tween the two wells. For the iV-atom occupation (filling factor y ) in the double-well trap 
we introduce the pseudo-angular momentum operators defined by 



PHASE TRANSITION 




(5) 




(6) 



and 



S z 



- (n L - n R ) 



(7) 



with the total angular momentum 




The Hamiltonian Eq.(2) can be written as 



H = -2JS X - AS Z + K^l + K 2 S 2 X , 



(8) 



with the parameters given by K\ — Uq — Ui, K 2 = 2U2, which is nothing but a LMG model 
Hamiltonian [llj]. Here we adopt a new method to study the energy spectrum and related 
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QPT of Hamiltonian Eq.(8). To this end we begin with the Schrodinger equation 

H$(<j>) = E<S>((f)), (9) 

where E is the energy eigenvalue to be determined and the generating function is constructed 
in terms of the spin function of S z representation such as 

m = ± C ™ e^t (10) 

m =-s J{s — my.{s + my 

in which the pseudo-spin operators have the form of differential operators 

S x = s cos (f) — sin0— , 
d(p 

S v — s sin (f) + cos 0-77; 

S --<A 

bz ~ l d<f>- 

Obviously S z and (p are two canonical variables describing the atom-number and phase 
differences respectively between the double wells. In the following we may convert our 
eigenvalue problem Eq.(9) to an effective single-particle Hamiltonian by making use of a 
proper unitary transformation and introducing an incomplete elliptic integral coordinate 
such that 

x= f* # =F(0,A), (11) 
Jo yd -A 2 sin 2 

where F{<p, A) is the incomplete elliptic integral of the first kind with modulus A 2 = K 2 /Ki, 
and then can achieve an effective single-particle in a quasi-periodic potential V(x) for the 
symmetric double- well A = 0, 



$(x) = Ef(x), (12) 



V(x) ^ a .^J^f, (13) 



with 

a(6 - 6) 2 - KCi - 6) + c 



v„ 



A' 2 + A 2 (6-6) 



2 



Lo A' 2 , l,aA' 2 
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and 

_ b 
fi ~ 2(a-A 2 Kni„)' 

where cn(x) and dn(x) denote the Jacobian elliptic functions of modulus A 2 and A' 2 = 1 — A 2 . 
Three parameters a, b, and c are related to the model parameters given by 

a = K 2 ^ + l)-J*/K 1 , 

b=J(N + l), 

and 

J 2 

The effective-Hamiltonian eigenvalue-problem Eq.(12) corresponds to a single-particle of 
mass m = \j1K\ in an effective potential V(x) plotted in Fig.3, which possesses two degen- 
erate minima located at 

x± = ±cn~ 1 {fi) mod[4/C(A 2 )], (14) 

respectively, where /C(A 2 ) is the complete elliptic integral of the first-kind with modulus A 2 . 
The two degenerate minima are separated by the central barrier located at x = with the 
barrier height given by 

Knax=(a-AV min )(l-| /U |) 2 , 

which can be controlled by the Josephson coupling constant J. Particularly when the Joseph- 
son coupling constant vanishes J = 0, we have /i — . the degenerate minima are located 
at x±(fi = 0) = ±/C(A 2 ) mod[2/C(A 2 )] and the potential V(x) becomes periodic. Notice the 
relation between the elliptic-integral coordinate x and the phase-angle the two degener- 
ate ground states correspond to the degenerate ±| mod(7r)-phase states of cold atoms in 
the double- well trap (Fig.4(a)). The QPT has obvious meaning in a real spin system, in 
which Si (i = x,y.z) denote the collective spin-operators. The degenerate ±| mod(7r)-phase 
states are nothing but the ferromagnetic phase of long range magnetic order (gapless) with 
two degenerate equilibrium-orientations of the magnetization. When the Josephson cou- 
pling constant J increases the height of center potential-barrier decreases (see Fig. 4) and 
the two- fold degeneracy of the ground state gradually lifts. At a critical point (Fig.4(c)) 



(l c = ±1 



(15) 
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the barriers located at x bl = mod[4/C(A 2 )] (corresponding to = 0mod(27r)) and x b2 = 
2/C(A 2 ) mod[4/C(A 2 )] (corresponding to = 7rmod(27r)) vanish, and we have the disordered 
phase called the paramagnetic phase with a energy gap in the spin language, which in our 
case of could atoms in a double-well trap indicates zero-and 7r-phase oscillations. The QPT 
from the long-range magnetic order to disordered phases is of the second-order coinciding 
with the Landau second-order phase transition theory, where the dimensionless barrier height 

'^(a-AH-r' 1 - 1 " 1 ' 
may be chosen as the order-parameter. The highest value of order-parameter h = 1 (/z = 0) 
corresponds to the degenerate ±| mod(7r)-phase states i.e. the ferromagnetic phase of two- 
fold degeneracy and the magnetic order decreases with the decreasing order-parameter h 
(increasing Josephson coupling constant J ). When the order-parameter vanishes h c = 
(|/x c | = 1) by properly adjusting the parameters J and N, the system approaches the 
non-degenerate mod (2n) -phase ( /i c — 1) or 7rmod(27r)-phase ( \i c = —1) ground state 
i.e. the paramagnetic phase in the spin language (Fig.4(c)). The phase-plane portraits 

2 

corresponding to the classical Hamiltonian of Eq.(12) H = ^ +V(x) are plotted in the 
Fig. 4 in comparison with the variation of potential barrier. The closed orbits show the 
oscillations with fixed average phase- difference between the double-well trap and waved- 
open lines indicate atom-number self-trapping. The variation of energy gap as a function 
of J/Uq with various K 2 is plotted in Fig.5 showing clearly the smooth second-order phase 
transition. It is worthwhile to remark that the entirely new method, which converts the 
many-body system to an effective single-particle in a potential field, has obvious advantage 
to study the dynamics and QPT of a complex system in a simple and visible way. 



IV. QUANTUM TUNNELING AND ATOM-NUMBER PARITY EFFECT 



In the previous section the dynamics of cold atoms trapped in a double-well potential 
has been investigated in terms of two conjugate variables, i.e. the atom-number-occupation 
difference and the relative phase in the elliptic integral coordinate x. In the case of large 
atom-number N quantum tunneling effect becomes negligibly small, however for small N, 
the quantum tunneling has to be taken into account. We demonstrate an atom number- 
parity effect of tunneling resulted from the boundary condition of the wave function $(0), 



12 



which is obviously periodic for even-number of N (and therefore the pseudo-spin s = y is 
integer) and antiperiodic for odd-number 

$(0 + 2tt) = $(0)e l2 ™, (16) 

while the boundary condition of wave function ip(x) in the elliptic integral coordinate x is 

^(x + 4/C) = (-l) 2s ^{x). (17) 

We begin with the transition amplitude, i. e. the Feynman propagator, between two degen- 
erate vacua x± induced by the quantum tunneling, which can be studied by the instanton 
technique 

< x+(T)\x-(-T) >= J V{x}e~ SE (18) 

where 

S E = f_ T C E dr (19) 

is Euclidian action defined with the imaginary time t — it and the Euclidian Lagrangian is 
given by 

C E = \m{^-f + V{x). (20) 

The Euclidian Feynman propagator can be evaluated in terms of stationary-phase 
perturbation-method, in which the zero-order perturbation comes from the action of clas- 
sical trajectory of pseudo-particles in the barrier region called the instantons. The explicit 
instanton solution of the classical equation of motion 

d 2 x _ dV(x) 
dr 2 dx ' 

is found in our case as 

, / 2r7tanh(^) \ 

xi(r) = tn- 1 '- , (21) 

\l- 77 2 tanh 2 (f ) J 

i / 2^ 1 tanh(f) \ 

x Jt) = tn' 1 ' , (22) 

1 ^ ^l-7T 2 tanh 2 (f)y 1 ^ 

which exist in the large (counterclockwise rotation) and small (clockwise rotation) barriers 

respectively (see Fig. 3), where tn{x) = sn(x)/cn(x) is the Jacobian elliptic function, r\ = 

^(1-^/(1 + /!), and 

uj = ^4tfi(a- A 2 V min ) (23) 
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is the oscillation frequency at the bottom of the potential well. Taking into account of 
contributions of both instantons X\{t) and x 2 (t) with the corresponding boundary conditions 
Eq.(17) the Feynman propagator is obtained up to the one-loop approximation and then 
the ground state tunnel splitting is abstracted as 



Ae = Qe' p \j2 (cosh 7 + cos 2ns) (24) 

where 

2 5 / 2 (a - AV min ) 3 / 4 



P = V ifiA 2 ln i-Av/l-^ ' 



and 



2/i /a - A 2 Km„ + A'v^T? 
7 = —7- W — arctan ■ 



A' V A': // 
To see the particle-number parity effect clearly we consider the case of vanishing Josephson 

coupling constant J = and thus /i = 7 = 0. For even-number of atoms (s is integer) the 

tunnel splitting is 

Ae = 2Qe~ p , 

however the tunnel splitting vanishes for odd-number of atoms (s is half-integer) called the 
quench of quantum tunneling due to the quantum phase interference of tunnel paths with 
the antiperiodic boundary condition of wave function Eq.(17). Thus the degeneracy of ±|- 
phase states cannot be removed by quantum tunneling. The tunnel splitting as a function 
of J is shown in Fig. 6 for odd- (a) and even- number (b) N respectively. 



V. CONCLUSION 



We conclude that the BH Hamiltonian in the superstrong interaction regime ought to 
be extended to include the two-body interaction of nearest neighbors, which results in a 
fundamental phenomenon of many-body system, namely the atom-pair tunneling. New dy- 
namics of various oscillations depending on the competition between the Josephson coupling 
constant J and interaction constant Uq is found in superstrong interaction regime, where the 
BH Hamiltonian gives rise to the fixed atom-occupation-number state only. The QPT and 
the critical transition point are analyzed analytically in terms of the potential-field method, 
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which allows us to convert the system of N cold-atoms in a double-well trap to an effec- 
tive single-particle in a quasi-periodic potential in the elliptic integral coordinate. The new 
oscillation states and related QPT should be observed in practical experiment. 
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APPENDIX A: 

The two-body interaction Hamiltonian in the second-quantization formulation is 

H in t = -J ^(X 2 )¥(X 1 )U(\ x 1 — x 2 \)^/(xi)^(x 2 )dx 1 dx 2 , (Al) 
with ty(x) being the field operator of ultracold atom-gas clouds, which in the lattice-mode 



expansion |l2l. Il3l | such that 

V(x) = J2 a i w i(%) (A2) 

i 

becomes 

1 Nl r 

Hm = - a \ a ) a k^i J w i (x 2 )w j (x 1 )U(\x 1 - x 2 \)w k (x 1 )wi(x 2 )dx 1 dx 2 , (A3) 

i,j,k,l 

where four sums over the lattice sites are independent and Nl denotes the total number of 
lattice sites. In the strongly interacting regime the sum over lattice sites should be extended 
to include nearest neighbors. The arbitrary sum over lattice sites can be obtained up to the 
nearest-neighbor approximation as 

Hint ~ Hint-o + Hint-n (A4) 

where the first term with all four modes on one-site 

Uo Nl ■ Uq Nl 

Hint-o = — Yl a l a l a i a i = yE"'^- !) ( A5 ) 

i i 

is the well known on-site approximation in the BH Hamiltonian. The sum over nearest- 
neighbor can be grouped under two configurations: 

Hint—n = Hint—nC^i 2) + Hi n t—niXi 
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where H int _ n (2,2) denotes any two-mode on one-site and H int _ n (l, 3) denotes three-mode 
on one-site. Obviously 

N L jj N l 

H int _ n (2, 2) = (C/i + U 2 ) J2 n i n i+i + ^2{.a\a\a i+1 a i+1 + a\ +l a\ +l aiai) (A6) 

and 

U 3 Nl 

Hint~n{^-, 3) = — ^(Oja i+1 aj + ia i+ i + a^a^i^a^i + a i+1 a i+1 ajaj + i + a i+1 a i+1 a i+1 aj 

i 

= U 3 J2( n i+i +rk- \){a\a i+ i + a| +1 Oi), (A7) 

i 

which can be combined with the hopping term in the Hamiltonian Eq.(l). 
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Figure Caption: 

Fig.l (color online)The time-evolution of average position for two-atom occupation in 
the weakly (J/U = 1.5) (a) and strongly (J/U = 0.2, 0.1) (b,c) interacting regime. Black 
dots denote the experimental data, red solid-line is the value evaluated from the proposed 
Hamiltonian and blue dash-line shows the result from Bose-Hubbard Hamiltonian. 

Fig. 2 (color online) Average position (a), visibility (b) and phase (c) for single-atom occu- 
pation in strong interaction regime (J/U = 0.2) with black dots denoting the experimental 
data. Coinciding red (solid) and blue (dash) lines are the values evaluated from the proposed 
and Bose-Hubbard Hamiltonians respectively. 

Fig. 3 Effective potential in elliptic-integral coordinate x (unit of /C) with asymmetric 
twin barriers and degenerate minima located at x± . x\ and X2 denote the clockwise and 
counterclockwise tunnel paths. 

Fig. 4 (color online) Phase-plane portraits for // = (a), 1 > /x > (b), and /i — 1 (c) and 
the variation of potential with respect to /i. 

Fig. 5 (color online) The energy gap in unit of K\ as a function of J evaluated by the 
numerical diagonalization of the Hamiltonian Eq.(8) with A = for the particle number 
N = 20 and various values of K 2 . 

Fig. 6 (color online) The energy gap as a function of J evaluated by the numerical diago- 
nalization of the Hamiltonian for the particle number N = 7 (a) and N = 8 (b). 
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